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We report on the experimental quantum teleportation of strongly 
nonclassical wave packets of light. To perform this full quan- 
tum operation while preserving and retrieving the fragile non- 
classicality of the input state, we have developed a broadband, 
zero-dispersion teleportation apparatus that works in conjunc- 
tion with time-resolved state preparation equipment. Our ap- 
proach brings within experimental reach a whole new set of hy- 
brid protocols involving discrete- and continuous-variable tech- 
niques in quantum information processing for optical sciences. 

In the early development of quantum information processing (QIP), a commu- 
nication protocol called quantum teleportation was discovered (/) that involves 
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the transportation of an unknown arbitrary quantum state \ip)by means of entan- 
glement and classical information. Experimental realizations of quantum telepor- 
tation (2, 3) and more advanced related operations (4) in the continuous-variable 
regime have been achieved by linear optics methods, although only for Gaussian 
states so far. However, at least third order nonlinear operations are necessary 
for building a universal quantum computer (5) — something that Gaussian oper- 
ations and Gaussian states alone cannot achieve. Photon subtraction techniques 
based on discrete-variable technology can provide useful nonlinearities and are 
used to generate Schrodinger's-cat states and other optical non-Gaussian states 
(6). Schrodinger's-cat states are of particular interest in this context, as they have 
been shown to be a useful resource for fault-tolerant QIP (7). It is therefore nec- 
essary to extend the continuous-variable technology to the technology used in the 
world of non-Gaussian states. 

We have combined these two sets of technologies, and here we demonstrate 
such Gaussian operations on nonclassical non-Gaussian states by achieving ex- 
perimental quantum teleportation of Schrodinger's-cat states of light. Using the 
photon- subtraction protocol, we generate quantum states closely approximating 
Schrodinger's-cat states in a manner similar to (8-11). To accommodate the re- 
quired time-resolving photon detection techniques and handle the wave-packet 
nature of these optical Schrodinger's-cat states, we have developed a hybrid tele- 
porter built with continuous-wave light yet able to directly operate in the time do- 
main. For this purpose we constructed a time-gated source of Einstein-Podolsky- 
Rosen (EPR) correlations as well as a classical channel with zero phase dispersion 
(12). We were able to bring all the experimental parameters up to the quantum 
regime, and we performed successful quantum teleportation in the sense that both 
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our input and output states are strongly nonclassical. 

A superposition of the quasi-classical coherent state \a) is one of the con- 
sensus definitions of a Schrodinger's-cat state |cat), typically written |cat) oc 
| a) ± | — a). Such optical Schrodinger's-cat states are known to be approx- 
imated by multiple photon subtractions from a squeezed vacuum state (6). In 
these protocols, a squeezed vacuum state S'(s) |0) is weakly tapped via a subtrac- 
tion channel, where |0) is the vacuum state and S(s) the squeezing operator with 
squeezing parameter s. When a photon detection event occurs in the subtraction 
channel, S(s) |0) is projected by the quantum action of the photon detector onto a 
non-Gaussian state, which can be tuned to approximate a Schrodinger's-cat state 
(8-10). The approximation is not perfect and can be quantified by means of the 
fidelity figure F cat = |(cat|a£(s)|0)] 2 (13). 

To represent the superposition nature of these states, we use the Wigner for- 
malism where for any quantum state |0) one associates a quasi-probability dis- 
tribution W(x,p), where x and p are the phase-space position and momentum 
parameters. W(x,p) is called the Wigner function and holds information exactly 
equivalent to \<f>) (14). Although the position x and momentum p quadratures 
operators and the vector state |0) are abstract objects, W(x,p) is always a defi- 
nite real- valued function that can be numerically reconstructed if one performs a 
complete phase-resolved sequence of homodyne measurement x cos 9 + p sin 9, 
a process called quantum tomography (15, 16). W(x,p) is not a true probability 
distribution, however, as there exist quantum states whose Wigner functions are 
not positive. \<j>) is defined to be a strongly nonclassical state when its Wigner 
function W(x,p) fails to be a positive distribution. Negativity in W(x,p) turns 
out to be an especially useful description of the nonclassicality of a Schrodinger's- 
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cat state |cat); \a) and | — a) induce two "classical" Gaussians in phase space, the 
superposition of which creates an oscillating interference pattern inducing nega- 
tivity in W(x,p). In contrast, a statistical mixture of \a) and | — a) would never 
show such negativity. 

In a quantum teleportation process, the input W m and output W ov x Wigner 
functions are related by the convolution (noted o) 

W 0Ui = W- m 0G e -r (1) 

where r is the EPR correlation parameter, G a is a normalized Gaussian of standard 
deviation a, and h (Plank's constant divided by 2ir) has been set to 1 (17). When 
finite quantum entanglement r is used, W ou i will be a thermalized copy of W m . 
Only with infinite r will G a become a delta function so that W m = W 0Vi [. The 
quality of quantum teleportation is usually evaluated according to the teleportation 
fidelity F tclc = (0 in |p out |0 in ), which can be written as F tele = 1/(1 + e~ 2r ) for 
Gaussian states (18). More importantly for our case, negative features of W m (if 
any) can only be teleported and retrieved in W OV £ when F te i e > 2/3 (19), a thresh- 
old also known as the no-cloning limit (20). However, the practical lower bound 
on Fteie will be higher because of decoherence and experimental imperfection of 
W m (21). We have thus defined the success criterion of Schrodinger's-cat-state 
teleportation as the successful transfer of its nonclassical features, or alternatively, 
successful teleportation of the Wigner function W m negativity. 

Our experimental quantum teleporter and Schrodinger's-cat-state source (Fig. 
1) upgrade the experiments described in (3) and (10), respectively. We use three 
optical parametric oscillators to generate the necessary squeezed vacua. One is 
used for the Schrodinger's-cat-state preparation; the other two are combined to- 
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gether on a half beam splitter whose two exit ports are the resulting pair of EPR- 
correlated beams. The teleportation is conducted in three steps. Alice first receives 
both the input state and her EPR beam and performs two joint quadrature mea- 
surements, obtaining results x and p . Bob then receives Alice's measurements 
(3 = (xo + ipo)/\/2 through the classical channels and applies the displacement 
operator D(f3) on his EPR beam. A final stage consists of a third homodyne detec- 
tor for tomography at the teleporter output. We emphasize that Alice and Bob do 
not assume any prior knowledge of the input state and adhere to unity-gain tele- 
portation, so that the teleporter does not have any restriction regarding the specific 
family of quantum states it can faithfully teleport. 

To benchmark our teleporter, we first evaluate the fidelity F te i e of teleporta- 
tion of the vacuum state |0), the coherent state of amplitude zero. At quantum 
optical frequencies where the mean thermal photon number is virtually 0, this 
is simply done by blocking the input port of the teleporter. The teleported vac- 
uum photocurrent is expected to have uniform Gaussian statistics with a variance 
a 2 = 1/2 + e~ 2r (h = 1) from which we can deduce teleportation fidelity (Fig. 2). 
The blue traces are the shot-noise level, the noise spectrum of the input vacuum 
|0). The red traces are the classical limit of teleportation obtained by turning off 
the entanglement between Alice and Bob (r = 0). We measure 4.8 dB of added 
noise above the shot noise, in agreement with the expected teleportation fidelity of 
0.5. When Alice and Bob share entanglement, the added noise drops to that shown 
by the green traces: 1.4 dB above the shot noise around 1 MHz, corresponding to 
a fidelity of 0.83. This is in agreement with the experimental figure of —6.9 dB 
that we observe in direct measurement of the EPR correlations shared between 
Alice and Bob. 
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In contrast to quantum teleportation experiments conducted to date for narrow 
sidebands of light (2, 3), our setup operates over a wide frequency bandwidth, 
as required by the nature of our input state. Because its generation relies on the 
detection of a single photon and the induced projection, a Schrodinger's-cat state 
made via photon subtraction is a short wave packet of light. A phenomenological 
way to picture these wave packets is to consider them as the closed boxes contain- 
ing the macroscopic superposition states as in Schrodinger's original idea. This 
requires Alice and Bob to teleport every frequency component of these "box-like" 
wave packets for faithful teleportation to occur. In this way, Alice and Bob do not 
need to actually teleport the Schrodinger's-cat states directly, but merely the po- 
tential boxes containing them. Consequently, Alice and Bob do not need to know 
when a detection event occurs; rather, they are only concerned with continuous 
and faithful "box" wave-packet teleportation, whichever state lies in the box. In 
fact, Alice and Bob actually teleport most of the time a squeezed vacuum state 
S(s)\0). 

In essence, our teleporter is a time-resolved apparatus that deconstructs the 
input wave packets into a stream of infinitely small time bins and reconstructs 
them at the output, within the extent of what we refer to as the teleportation band- 
width. This bandwidth is clearly visible in both green experimental traces where 
the added noise slowly increases with frequency (Fig. 2). This is a direct con- 
sequence of the finite bandwidth of squeezing used for entanglement. However, 
across the frequencies relevant to our input state, teleportation fidelity is always 
greater than the no-cloning limit of 2/3, a necessary regime for negativity telepor- 
tation. A very careful implementation of the classical channel has been required 
(12) to achieve experimental realization of this fidelity. 
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To verify the success of Schrodinger's-cat-state teleportation, we perform ex- 
perimental quantum tomography of the input and output states independently (Fig. 
3). Both input and output marginal distributions exhibit the characteristic eye 
shape of photon- subtracted squeezed states, with a clear lack of detection events 
around the origin for any phase. Although necessary, this feature alone is not suf- 
ficient to confirm the presence of negativity in W m or W OU [. The reconstructed 
input Wigner function W^ n shows the two positive Gaussians of \ct) and | — a) 
together with a central negative dip (W m (0, 0) = —0.171 ± 0.003) caused by the 
interferences of the |a) and | — a) superposition. The output Wigner function 
W ou t retains the characteristic non-Gaussian shape as well as the negative dip 
(W out (0, 0) = -0.022 ± 0.003) to a lesser degree. The degradation of the central 
negative dip and the full evolution of W m towards W OU [ can be fully understood 
using Eq. 1 with a model of W m , as was done in (21). Given the measured input 
state negativity of W^O, 0) = —0.171 and —6.9 dB of squeezing, the results of 
(21) predict an output negativity value of —0.02, in good agreement with mea- 
sured output negativity. Although this figure does not take into account the input- 
state squeezing, a more detailed model shows that a squeezing parameter s = 0.28 
affects output negativity in the third decimal place only (12). The experimental in- 
put and output states have an average photon number (n) equal to 1.22 ± 0.01 and 
1.33 ± 0.01, respectively (12). The increase in the output-state size is due to tele- 
portation induced thermalization. We calculate that the fidelity F cat is as high as 
0.750 ± 0.005 for the input Wigner function W m , with the nearest Schrodinger's- 
cat having an amplitude (ami 2 = 0.98 (12). However, after the teleportation W OU [, 
fidelity is reduced to 0.46 ± 0.01, with the nearest Schrodinger's-cat state having 
an amplitude |« ou t| 2 = 0.66. If W ou i fidelity is calculated with |a in | 2 = 0.98, 
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thenF cat = 0.45 ±0.01. 

We have demonstrated an experimental quantum teleporter able to teleport full 
wave packets of light up to a bandwidth of 10 MHz while at the same time pre- 
serving the quantum characteristic of strongly nonclassical superposition states, 
manifested in the negativity of the Wigner function. Although F cat and W(0, 0) 
drop in the teleportation process, there is no theoretical limitation other than avail- 
able squeezing, and stronger EPR correlations would achieve better fidelity and 
negativity transmission. The various more complex states generated as an applica- 
tion of photon- subtraction so far (22, 23) can be, in principle, readily sent through 
our broadband quantum teleporter. This opens the door to universal QIP and fur- 
ther hybridization schemes between discrete- and continuous-variable techniques 
(24). 
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Figure Captions 
Figure 1 

Experimental setup. OPO: Optical Parametric Oscillator, APD: Avalanche Photo 
Diode, HD: Homodyne Detector, LO: Local Oscillator, EOM: Electro-Optical 
Modulator, ADC: Analog to Digital Converter, FC: Filtering Cavity. See (12) for 
details. 

Figure 2 

Broadband teleportation of the vacuum state |0). Experimentally measured 
power spectra of the measured photocurrents calculated by Fourier transform are 
shown for the position (A) and momentum (B) quadratures, (blue) shotnoise in- 
put, (green) quantum teleportation output, (red) teleportation output without en- 
tanglement. Reconstructed Wigner functions of the input state |0) (C), quantum 
teleported vacuum (D) and classically teleported vacuum (E). 

Figure 3 

Teleportation of Schrodinger's cat states. Experimentally measured input state's 
Wigner function W m (A), marginal distribution (B) and photon number distribu- 
tion (C), as well as output state's Wigner function W ov x (D), marginal distribution 
(E) and photon number distribution (F). 
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Classical Channels 
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Figure 3: 
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Supporting Online Material 

1 Experimental setup 

We use a continuous-wave titanium- sapphire laser at 860 nm (SolsTiS-SRX, M 
Squared Lasers, Verdi 1 8 W pump, Coherent) delivering a round, zero-order Gaus- 
sian TEM beam. 1 W of that beam is used to pump a second-harmonic generation 
cavity with a KNb0 3 crystal for squeezing generation. Weak coherent beams 
are injected into each OPO and are phase controlled with respect to each other 
by monitoring classical phase interference signals. To support phase control and 
noiseless measurement at the same time and allow for time-gated operations, these 
beams are alternatively switched on and off using acousto-optical modulators at 
a frequency of 2 kHz. Phase locks are achieved when the beams are switched on 
while measurement is performed when the beams are switched off. 

2 Schrodinger's cat states source 

We use for the photon subtraction protocol a weakly squeezed vacuum S\0) gener- 
ated with an Optical Parametric Oscillator (bow-tie cavity, PPKTP crystal, finesse 
48, Full-Width-at-Half-Maximum 12.3 MHz). We pump this OPO with 15 mW 
of 430 nm light, for a measured squeezing parameter s = 0.28 ± 0.02. Intra- 
cavity losses with the pump injected are 0.0033 ± 0.0005. From the produced 
squeezed vacuum a small fraction (3%) is tapped on a beam-splitter and guided 
into an Avalanche Photo Diode (APD, Perkin-Elmer). Two Fabry-Perot cavities 
(FWHM: 1 10 and 36 MHz, finesse: 1540 and 1630) on the way to the APD allow 
to select the degenerate central mode of OPOl only. The photon detection rate is 
9000 ± 500 events per second. The dark count rate is 120 ± 20 events per sec- 
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ond, for a event-to-dark-counts ratio of 66 ± 6. The actual values of R and s are 
the result of optimization between event rate and squeezing purity. Rather than 
maximizing the overlap with the nearest cat state, we are actually more concerned 
in the input state negativity W^O, 0) since our criterion of success is negativity 
Wo U t(0, 0) of the output teleported state. By lowering the pump beam power, the 
purity of the squeezing is increased which is beneficial for the input state cen- 
tral negativity W^ n (0, 0). Furthermore it reduces the probability of higher order 
pairs of photons which is also beneficial for W^O, 0). However it decreases the 
event-to-dark-counts ratio, which increases W-^O, 0). 

3 Broadband teleportation 

Compared to the bandwidth of OPOl used for the generation of Schrodinger's cat 
states, we use broader bandwidth OPOs for the EPR entanglement source. Both 
OP02 and 3 are built on bow-tie cavities with PPKTP crystals. FWHM is 24 MHz 
and finesse is 41. Pump beam powers are 120 mW and intra-cavity losses with 
pump beam injected are respectively 0.0035 and 0.0037. 

For teleportation to work on a broadband range of frequencies, every elec- 
tronic stage in the classical channel needs a flat gain band as well as a linear phase 
delay. Only in that case will classical information received by Bob from Alice 
match in time and in magnitude the EPR correlations shared between them at ev- 
ery frequency. A nonlinear phase delay across the teleportation bandwidth would 
distort the wave-packet shape while a non-flat gain would result in excess noise at 
specific frequencies. 

Bob performs displacements by mixing two auxiliary coherent beams with 
his EPR correlated beam on two 99.5/0.5 beam splitters to preserve quantum 

15 



correlations. Each auxiliary beam is phase modulated according to Alice's results 
for the x and p quadratures. To accommodate with the DC- 10 MHz bandwidth of 
teleportation, we use for modulation of these auxiliary beams optical waveguide 
phase modulators based on Mach-Zehnder interferometers with a matched 50 
impedance that allows for easy broadband operations. 

To compensate for the global time delay of the slower electrical classical chan- 
nels, the optical path length of Bob's EPR beam is set at 12 m which allows syn- 
chronized overlap between classical information and EPR correlations owned by 
Bob at the output of the teleporter, both detrimental to the quality of teleportation. 

Finally the gain of the classical channels is adjusted to unity following the 
method of (3). 

4 Homodyne detection and tomography 

The three homodyne measurements are performed in the usual way with 50/50 
beam splitters. For better spatial mode-matching the local oscillator beams are 
filtered in a separate bow-tie cavity. Visibilities are simultaneously above 99% for 
every pair of signal beam and LO beam. Homodyne detectors are built according 
to the transimpedance design and employ a pair of silicon pin photodiodes with 
a 99% quantum efficiency (Hamamatsu, special order). The bandwidth of the 
detectors is above 20 MHz. Local oscillator power is set at 10 mW for every 
homodyne measurement to avoid saturation of the photodiodes yet to obtain a 
high signal-to-noise ratio. 

For quantum tomography, the photocurrent to be analyzed is sent to an analog 
to digital converter (National Instruments, PXI-5152, sampling rate 1 GHz, 8 bits 
resolution). To retrieve the non-Gaussian statistics of the state the AD converter 
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is triggered by the clicks of the APD and records 500 data points for every trigger. 
These 500 points together constitute a single quadrature measurement calculated 
using the temporal mode function of our input state taken as exp [— j\t\) where 7 is 
the bandwidth of OPOl. One set of measurements for a full tomography is made 
of 200000 quadrature measurements with a uniformly distributed phase angle 9 
and recorded in approximately 90 seconds. All Wigner functions are reconstructed 
using the maximum-likelihood iteration algorithm without any compensation for 
imperfection or prior knowledge of the quantum state. 



5 Estimation of photon size and Schrodinger's cat 
state fidelity 



From the experimentally reconstructed density matrices, we calculate the input 
and output state average photon number (n) using the formula 

N 

(n) = tr (rip) = n Pnn (3) 

n=l 

where N + 1 is the photon size of the Hilbert space used in the maximum- 

J1 

likelihood iteration algorithm (N = 1 5LJ). To calculate the Schrodinger's cat state 
fidelity F cat = (cat|p|cat) we express the scalar product in phase space as 



where is the Wigner function associated to the density matrix p. Our target 
Schrodinger's cat state is the so-called odd-cat state (\a) — | — a))/ y/N where N 
is a normalization factor expressed as 2(1 — exp [—2a 2 ]) and a is a real number. 

'For N above 10 there are virtually no differences in the reconstructed states and the numerical 
differences have been checked to be much smaller than the experimental error boundary. 




(4) 
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Its Wigner function W ca , t is given as 

W cat (x,p) = ^e~ x2 ~ p2 (f 2al cosh(2v / 2ax) - cos(2>/2ap)) (5) 

To estimate the fidelity to the nearest odd-cat state, using our reconstructed density 
matrix p and the above expression of W cs x we evaluate numerically the integral 
(3) and look for its maximum value over the free parameter a. 

6 Evolution of negativity in the teleportation 

Using equation (1) we want to predict the expected negativity at the output of the 
teleporter given the EPR correlation parameter r and the input Wigner function 
W m . As W m is neither a pure state nor a Gaussian state there is no easy way to 
characterize it and we need to use a theoretical model to approximate the exper- 
imentally measured W m . The model proposed in (21) uses a mixture of vacuum 
|0) and one photon state |1) written 

Pl =7 7 |l)(l| + (1-7 7 )|0)(0| (6) 

If W^(x,p) is the Wigner function of p x , then the negativity of this state is ex- 
pressed by W^(0, 0) = (1 — 2r])/n. Using equation (1) we derive the negativity 
of the output state Wigner function, had p x been used as the input state: 

W£2 (0, 0) = (l - 277 + 2e- 2r ) /n (l + 2e~ 2r ) (7) 

Using the measured input state negativity we calculate r\ = 0.77. From the mea- 
sured squeezing of OP02 and 3 we deduce r = 0.795, from which we calculate 
w£>(0,0) = -0.0207. We need a more advanced model to include the effect 
of the squeezing of the input state (quantified by the squeezing parameter s) on 
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the output state negativity. For that we first think of the input state as the state 
\4>2) = aS(s)\ty whose Wigner function is written 

W i2) (x, p) = 2(e 2s x 2 + e' 2s p 2 - 1/2) exp [-e 2s x 2 - e~ 2s p 2 ] (8) 

However |0 2 ) is a pure state and does not model the experimental measured 
Wigner function W m well. Especially it has maximum negativity —I /it. We 
model W m using W^^x^p) together with the linear beam splitter loss model 
with loss parameter 1 — rj: 

H K»( I ,p) = I( W ».G»)(i,.2L) (9) 



with A = y ^ and G a a normalized Gaussian of standard deviation a. Using 
again equation (1) on we finally derive the output negativity 14^(0, 0) ex- 
pressed as 

lOO-O) - vru; '' -" -372 (10) 



7 r(^ 2 + 477(^-77)sh 2 (s)) 1 



with g r — 1 + 2e 2r . Using the same experimental figures as before and adding 
s = 0.28, we deduce r] = 0.79, r = 0.795 and W^ t (0, 0) = -0.0247. We notice 
that the output negativity is increased when using the second model of input state. 
This can be understood when looking at the higher value of r\ needed to fit the 
input negativity W^ n (0, 0). Therefore is actually a purer state than and 
is easier to teleport in term of negativity. 
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